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ABSTRACT 



Aims. The influence of a toroidal magnetic field on the dynamics of Rossby waves in a thin layer of ideal conductive fluid on a rotating 
sphere is studied in the "shallow water" magnetohydrodynamic approximation for the first time. 

Methods. Dispersion relations for magnetic Rossby waves are derived analytically in Cartesian and spherical coordinates. 
Results. It is shown that the magnetic field causes the splitting of low order (long wavelength) Rossby waves into two different 
modes, here denoted fast and slow magnetic Rossby waves. The high frequency mode (the fast magnetic Rossby mode) corresponds 
to an ordinary hydrodynamic Rossby wave slightly modified by the magnetic field, while the low frequency mode (the slow magnetic 
Rossby mode) has new and interesting properties since its frequency is significantly smaller than that of the same harmonics of pure 
Rossby and Alfven waves. 

Key words, magnetohydrodynamics (MHD) — waves 



1. Introduction 

The large-scale dynamics of planetary atmospheres is mostly de- 
termined by Rossby waves. These waves arise because of the lat- 
itudinal variation of the Coriolis parameter and are widely used 
in the geophysical context (Pedlosky 1987, Gill 1982). Rossby 
waves also can be of importance in solar and stellar astrophysics, 
particularly in a thin layer called tachocline that is believed to 
exist below the convection zone of solar-like stars (Spiegel & 
Zahn T992l G ough & Mclntyre [19981 Garaud|2002l Callv 120031 
Miesch |2005l ). The thickness of the tachocline is very small com- 
pared to the stellar radius and, therefore, the ordinary shallow 
water approximation can be easily applied, but the hydrody- 
namic (HD) Rossby wave theory needs to be modified in the 
presence of a large-scale horizontal magnetic field. The influ- 
ence of the horizontal magnetic field on the large-scale fluid dy- 
namics has been studied in the context of the Earth's liquid core 
using the two dimensional /3-plane approximation in Cartesian 
coordinates by Hide d 1966b . However, to study the plasma dy- 
namics over spatial scales comparable to the stellar radius re- 
quires to consider spherical coordinates. Magnetohydrodynamic 
(MHD) "shallow water" equations for the solar tachocline have 
been recently proposed by Gilman d2000| l and the dynamics of 
various "shallow water" MHD waves in the solar tachocline have 
been studied by Schecter et al. (120011 ) (see also De Sterck l200TT) . 
Large-scale Rossby-like waves are absent from their considera- 
tion as the /-plane approximation has been used. However, we 
should mention the recent work by Leprovost & Kim (120071) . 
which studies the influence of shear, Rossby, and Alfven waves 
on the transport properties of MHD turbulence on a /?-plane in 
the solar tachocline. 



Here we use the MHD "shallow water" equations in order to 
study the influence of a toroidal magnetic field on the dynamics 
of Rossby waves in a rotating sphere. First we use Cartesian co- 
ordinates and derive the dispersion relation of magnetic Rossby 
waves in the /?-plane approximation. Next, we solve the problem 
in spherical geometry, thus deriving the propagation properties 
of magnetic Rossby waves with a wavelength comparable to the 
radius of the sphere. 



2. Basic considerations 

Let us consider a thin layer of ideal conductive fluid on a sphere 
rotating with angular velocity Jin. The layer is taken as an in- 
compressible fluid with a rigid base and a free upper surface and 
is permeated by a horizontal uniform magnetic field. The un- 
perturbed uniform thickness of the layer, Hq, is smaller than the 
stratification scale height and so the medium density can be con- 
sidered uniform. This system differs from the classical shallow 
water system only by the presence of the magnetic field. Then, 
the "shallow water" MHD equations in an inertial frame can be 
written as (Gilman 2000 Schecter et al. 120011) 



d,B + (V-V)B = (B-V)V, 



d,\ + (V-V)V = (B-V)B - gVH, 

\np 

d t H + V-(.ffV) = 0, 



(1) 
(2) 
(3) 
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where V and B are the horizontal velocity and magnetic field, 
H is the thickness of the layer, p is the fluid density, V is the 
horizontal gradient and g is the gravitational acceleration. The 
divergence-free condition for the magnetic field, which arises 
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from the requirement that B is parallel to the upper free surface, 
takes the form (Gilman 120001 



V-(flB) = 0. 



(4) 



A general feature of the large-scale dynamics in such a sys- 
tem is that it does not significantly depend on the chosen geome- 
try (fiat, spherical or cylindrical) for equatorially trapped waves 
(Longuet-Higgins 1965 Pedlosky |1987l l. However, the consid- 
eration of spherical geometry is desirable for waves with wave- 
length comparable to the size of the sphere. Therefore, we first 
study the problem in the simpler Cartesian coordinates and then 
turn to the more complicate spherical geometry. 



3. Cartesian coordinates 

We consider a local Cartesian coordinate frame (x, y, z) in which 
the x axis is directed towards the rotation, the y axis is directed 
towards the north pole of the sphere and the z axis is directed 
vertically. 

Let us next consider that the unperturbed magnetic field, 
(B x , 0, 0), is directed along the x axis. Then, after linearizing 
Eqs. (l)-(3) their components are written in the rotating frame 

as 



du x 
~d7 

dUy 

~dt 

db x 

~dT 



dh 

~dl 



+ H 



B x db x dh 


(5) 


4np dx ^ dx' 


B x dby dh 


(6) 


4-np dx ^ dy' 


dby duy 
"' ~dl ~ Bx ~dx , 


(7) 




(8) 



where u x , u y , b x and b x are the velocity and magnetic field pertur- 
bations, h = H— Ho is the perturbation of the layer thickness and 
/ = 20o sin is the Coriolis parameter (with the latitude). For 
zero magnetic field this system transforms into the HD shallow 
water equations (Pedlosky |1987l l. 

Differentiation with respect to time of Eqs. (5)-(6) and using 
Eqs. (7)-(8) gives 



d 2 u x 

~df~ 

d 2 Uy 

~~dW 



dUy 

f ~d~T 



,dh 



c du x 

~dT 



dx 2 

d 2 Uy 

A ~dx 2 



+ C 2 



d 2 u 



d 2 Uy 



dx 2 dxdy 



d 2 u x d 2 u. 



( 

\ dxdy dy 2 



(9) 



(10) 



where va = B x j yjAirp and Co = ygHo are the Alfven and sur- 
face gravity speeds, respectively. 

We now perform a Fourier analysis of the form exp(-icot + 
ik x x) and after some algebra obtain 



d 2 u y 

~d~f 



c 2 



k 2 v 2 

c 2 



2 f 2 



k x co df 



C 2 (co 2 - k 2 x v 2 A ) (co 2 - k 2 x v 2 A ) dy 



inclusion of the magnetic field leads to the modification of the 
wave modes. 

At a given latitude, 0o, one can perform a Taylor expansion 
of the Coriolis parameter and retain the lowest order latitudinal 
variation off, which leads to (Pedlosky |T987l Gill [T9"82j l 



f = fo+py, 

where the parameter 

2Q.Q 
B = — — cos ©o 



(12) 



(13) 



(with Rq the radius of the sphere) plays a major role in the so 
called /?-plane approximation. Away from the equator By <s /o 
and therefore from Eq. (1 1) we readily get 



d 2 i, 



dy 2 



k 2 v 2 

c 2 



/ ,2 f2 

u Jo 



k x coB 



Cl(u 2 



(cj 2 -k 2 x v 2 A ) 



= 0(14) 



Thus, we can now perform a Fourier analysis of the form 
exp{ik y y), which gives the dispersion relation 



u A -[2k 2 x v 2 A+ f 2 



C'iK + fyW-C'kxBu 
k 2 A[k 2 x v 2 A + Cl(k 2 x + k 2 )] 



0. 



(15) 



This dispersion relation contains high and low frequency 
branches, which respectively correspond to magneto-gravity 
waves and to Alfven and Rossby waves. Note that for B — this 
dispersion relation transforms into that of /-plane MHD "shal- 
low water" waves (Schecter et al. 1200 11 1. 

We next concentrate in the case of small Alfven speed, i.e. 
va <k Co, such as corresponds to the interiors of solar-like stars. 
Then, the high frequency branch of Eq. (15) contains Poincare 
waves, whose dispersion relation is (Pedlosky |1987l l 

<u 2 =f 2 +C 2 (k 2 x + k 2 ), (16) 

while the low frequency branch yields the dispersion relation 

k x B 



2 

CO + 



kl + k 2 



■rCO 



■k 2 v 2 



0. 



(17) 



Note that the dispersion relation (17) was first obtained by Hide 
(119661 1 in the two dimensional case (see also Acheson & Hide 
I1973I I. This formula reveals some interesting properties. For 
short wavelengths, i.e. large k x , the last term in Eq. (17) dom- 
inates over the second one, which leads to the solution 



co = +k x v A . 



(18) 



This is the dispersion relation of pure Alfven waves unaf- 
fected by rotation and propagating eastward and westward in the 
toroidal direction. 

Nevertheless, for large-scale motions pure Alfven waves no 
longer exist and instead we have Rossby waves modified by the 
magnetic field. For large wavelengths, i.e. small k x , Eq. (17) has 
two different solutions. For the high frequency solution one can 
easily recover the dispersion relation of HD Rossby waves, 



= ©11) 



k x B 
kl+kl' 



(19) 



When va — this equation governs the linear dynamics of vari- 
ous kinds of waves (namely Poincare, Kelvin and Rossby waves) 
in the HD shallow water approximation (Pedlosky 1987), but the 



For the low frequency solution we have the dispersion relation 

(20) 



k x v 2 A (kl + k 2 ) 
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Hence, the horizontal magnetic field causes the splitting of 
ordinary large-scale Rossby waves into two modes propagat- 
ing in opposite directions. The high frequency mode has the 
properties of HD Rossby waves and can be called fast mag- 
netic Rossby mode. But, additionally, a lower frequency mode 
arises whose frequency is significantly smaller than that of pure 
Alfven and Rossby waves at the same spatial scale. Due to its 
small frequency it can be called slow magnetic Rossby mode 
("hydromagnetic-planetary waves" in Acheson & Hide I1973K 
The phase speed of the mode in the x direction depends on both 
the Alfven speed and the B parameter, 



co v 2 A (k 2 x +k 2 y ) 



V p h = — = 



P 



(21) 



The phase speed is different from Alfven and Rossby phase 
speeds, which again indicates the different nature of this wave 
mode. 

Numerical dispersion diagrams for the general dispersion re- 
lation (15), i.e. without assuming va <k Co, are presented in 
Fig. 1 . The upper panel displays all wave solutions and shows 
that Poincare waves are almost not affected by the magnetic 
field. The middle panel is a detailed view of the low frequency 
branch of the dispersion diagram. It is clearly seen that for small 
scales, i.e. for large k x , magnetic Rossby waves tend to the 
Alfven wave solutions (dashed lines), whereas for small k x , i.e. 
for large spatial scales, the two modes behave differently: the so- 
lution with higher negative frequency corresponds to HD Rossby 
waves (triangles) and the low frequency solution, which differs 
from pure Rossby and Alfven wave dispersion curves, is a new 
wave mode. Finally, the bottom panel shows the perfect fit be- 
tween the solutions to Eq. (15) and the approximate dispersion 
relation for magnetic Rossby waves (circles), Eqs. (19)-(20), in 
the limit of small k x . 

The present consideration in Cartesian coordinates gives the 
basic properties of magnetic Rossby waves. Nevertheless, to 
study the dynamics of Rossby waves with spatial scales com- 
parable to the radius of the sphere it is desirable to use spherical 
coordinates. Therefore, in the next section we study the same 
problem in spherical coordinates (r, 9, <f>). 

4. Spherical coordinates 

Let us consider an unperturbed toroidal magnetic field Ba. Then, 
the linearized form of Eqs. (l)-(3) can be rewritten in the rotat- 
ing frame as 



5,4 



du e or> g dh 

— 2£2 () cos 6ua + — 

dt * Rode 4npR sme 

_ Ba cos e , 



db„ 



A-npRo sin# 



(22) 



— — + 2Q cos eu e + — 

dt Rosinedcp 4npR de 



DA 



db$ 



Ba cos e . 



AnpRosine deb AnpRo sin 6* 



b e = 0, 



Rq sin e dcf> 



dh H d H ditA 

— + — (sin eu e ) + — = 0, 

dt Rosinede ' " 

db e _ Ba due _ 
dt Rq sin d<p 



(23) 
(24) 
(25) 
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Fig. 1. (a) Dispersion diagram of "shallow water" waves in the 
presence of a horizontal magnetic field. The two extreme up- 
per and lower solutions correspond to Poincare waves, which 
are almost unaffected by the magnetic field, whereas the two 
low frequency modes are magnetic Rossby waves, (b) Low fre- 
quency branch of the dispersion diagram, (c) Low frequency, 
large wavelength part of the dispersion diagram. In all panels 
solid lines correspond to the solutions to Eq. (15), triangles to 
the HD Poincare and Rossby waves (i.e. in the absence of mag- 
netic field), dashed lines to pure Alfven waves (whose dispersion 
relation is given by Eq. (18)) and circles to the approximate an- 
alytical formulas (19)-(20). The parameters used to obtain the 
dispersion diagram are k y Ro - 1, Va/Cq = 0.005 and = 45°. 



dbA Id, v 
-7T- + —t-(u8Ba) = 0, 



Ro d6 



(26) 



where ug, ua, be and b,p are the velocity and magnetic field per- 
turbations, while h = H - Ho is the perturbation of the layer 
thickness. For zero magnetic field this system transforms into 
the HD shallow water equations (Longuet-Higgins 1965). 

We assume the unperturbed magnetic field to be Ba = 
Bo sin e, which means that it has a maximal value at the equa- 
tor and tends to zero at the poles. We take a sinusoidal depen- 
dence of the magnetic field on e for two main reasons: first, the 
sinusoidal profile simplifies the calculation in the spherical sym- 



metry and second, the toroidal magnetic field seems to be lo- 
cated mainly in low latitudes due to the eruption of magnetic 
flux at these latitudes. Let us next introduce the new variables 
tig = sinOug, u^ = sinOu^, bg = sindbg, b$ = sinOb^,. Then, 
Eqs. (22)-(26) take the form: 

dug ori , 8 . a dh B dbg 

— 212(1 cos Ouj, h sin 6— h 

dt R 80 4jrpR Q deb 
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and Eq. (38) takes the form 



fin 



du$ 
~dt~ 



+2Q.Q cos 9u,g + - 



dh 



AnpRo 
B db 



cos 6b$ = 
2B 



(27) 



• 2a dh H ° ■ n 

sin 8 — + — smO 

dt R Q 89 



R deb AnpR {) dtp 4npR 
dug Hq du^ 



Ro d(f> 



dbg 

dt 
db* 



go dug 

Ro~~tty 



Bo . dug 
— + — sin#— - 
dt Rq d9 



cos 6b e = 0, (28) 
(29) 
(30) 
(31) 



We now perform a Fourier analysis of the form exp(-ia>t 
iscb) and define 



2Clo~ 



A. 



gh 



4D. 2 R 2 



20()/?o 



cos ( 



- sin#- 



(l-/^ 2 )- 



d 



d6 v n 'dp 
After some algebra we get 
— A 2 Ug — pAu^ — ADrj + s 2 a 2 u g + 2a 2 pDUg = 0. 

0. 



D. 



A Uj, 



■ pAiig + Asrj — 2sa 2 pug — a 2 sDug 



6/1(1 — p )r]— Dug - 5M0 = 0, 



(32) 

(33) 
(34) 
(35) 



where iig = iug. 

Substitution of u$ from Eq. (34) into Eqs. (33) and (35) leads 

to 



9 

A ug + p u^ - {AD + sp)rj + s a Hg + 2a pDiig + 



ps—(D + 2p)Ug = 0, 
A 



eA(l -p )T]~ Dug + -(pug 
A 



«7) + ^2 r(D + 2fi)u e = 0. 



(36) 
(37) 



After obtaining rj from Eq. (37) and substituting it into Eq. (36) 
we get a single equation for u g , 



(AD + sp) 



1 



s 2 -eA 2 (l -p 2 ) 



AD- sp - - I s 2 A{D + 2p) 



OO OO O 

-(A — p )ug + s a Ug + 2a pDiig + ps — (D + 2p)ug 

A 



tin 



Q38) 



In the approximation for slowly rotating stars (such as in the 
solar case), 



e 4Q 2 R 2 



gH s 2 



« 1 



(39) 



[(AD + sp)(AD - sp) s 2 (AD + sp)(D + 2p) ■ 



s 2 (A 2 



■ p z ) + s 4 a 2 + 2a L s L pD + ps 3 —(D + 2p)]u e = 0. (40) 

A 



..2 „2. 



The approximation (39) implies that magneto gravity waves (i.e. 
magnetic Poincare and Kelvin waves) are neglected from our 
consideration, thus retaining only magnetic Rossby waves. 
Now, Eq. (40) can be rewritten as 



dp 



if 

n(n +1) 



dp 1 - p 2 

sA + 2s 2 a 2 
A 2 — a 2 s 2 



+ n(n + 1) 



llg = 0, 



(41) 



(42) 



Equation (41) is the associated Legendre differential equation 
(Abramowitz & Stegun 1964), whose typical solutions are the 
associated Legendre polynomials, 



u e = P s n (cos6) 



(43) 



if n is an integer (when n is not integer the solutions are the 
associated Legendre functions). 

Equation (42) defines the dispersion relation for spherical 
magnetic Rossby waves, 



(A\ 2 A 
-J + - +a 2 [2-n(n+ 1)] =0, 



(44) 



where n plays the role of the poloidal wavenumber. 

In the non-magnetic case, i.e. for a — 0, the dispersion rela- 
tion reduces to the HD Rossby mode solution (Longuet-Higgins 
11965 ). But the magnetic field causes the splitting of the ordinary 
HD mode into fast and slow magnetic Rossby modes as in the 
rectangular case (it is worth recalling that Rossby-type waves 
are sometimes called planetary waves, which is probably a more 
appropriate name for the spherical geometry). 

Before studying the dispersion diagram, let us first consider 
the properties of wave modes for particular values of n. For 
purely toroidal propagation, i.e. when n — 0, we have only one 
mode with the dispersion relation 



S 



(45) 



In the weak magnetic field limit this mode has significantly 
lower frequency than pure Rossby and Alfven waves and can be 
associated to the slow magnetic Rossby mode. Note, that n — 
mode is absent in non-magnetic spherical case (Longuet-Higgins 

rr%5i 

For n — 1 we have only the fast magnetic Rossby mode, 
which in this case is identical to the HD Rossby mode, with the 
dispersion relation (Longuet-Higgins l 1 9651 l 



CD 
S 



-Qo. 



(46) 



For n > 1 we have both fast and slow magnetic Rossby 
modes. In the weak magnetic field limit, i.e. for va <k 2Qnfio, 
the dispersion relation for the fast magnetic Rossby mode is 



CO 

s 



2Qo 
n(n + 1) 



(47) 
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Fig. 2. Numerical dispersion diagram of spherical "shallow wa- 
ter" waves in the presence of a toroidal magnetic field. The de- 
pendence of the wave frequency on the poloidal wavenumber n 
is plotted here for a 2 = 0.036. The continuous and dashed lines 
are the solutions for the slow and fast magnetic Rossby modes, 
respectively. The dotted line is the HD Rossby mode, which is 
obtained from Eq. (44) with a 1 = 0. 



Fig. 3. Frequencies of the s — 1, n = 2 harmonics of fast (dashed) 
and slow (continuous) magnetic Rossby modes vs. the ratio of 
the Alfven speed to the rotation rate, a = va/2£IoRq. When a 
tends to zero, then the slow magnetic Rossby mode has a very 
low frequency, while the fast mode frequency tends to ~ 0.3£2o- 



(which is similar to HD Rossby mode, Longuet-Higgins 1965) 
and the dispersion relation for the slow magnetic Rossby mode 



is 



[2-n(n + 1)]. 



(48) 



Equations (47)-(48) show that the toroidal magnetic field causes 
the splitting of ordinary HD Rossby waves into two different 
modes propagating in opposite directions for n > 1 . Such split- 
ting is also present for all non-integer values of n smaller than 
1 . Thus, the general behaviour of the wave modes in spherical 
geometry is similar to that of the Cartesian case. There are some 
differences, however. First, the dispersion relation of magnetic 
Rossby waves in a Cartesian frame depends on the latitude; sec- 
ond, magnetic Rossby waves are dispersive with respect to the 
toroidal wavenumber, k x , in Cartesian geometry (see Eqs. (19)- 
(20)), while they are not dispersive with respect to the toroidal 
wavenumber, s, in spherical geometry. 

Numerical solutions to the general dispersion relation (44) 
are presented in Fig. 2. The dependence of the wave frequency 



on the poloidal wavenumber n is plotted here for a 1 - 0.036. 
The continuous and dashed lines are the solutions for slow and 
fast magnetic Rossby modes, respectively. The dotted line is 
the HD Rossby mode, which is obtained from Eq. (44) with 
a 2 = 0. It is clearly seen that for small scales, i.e. for large n, 
fast and slow magnetic Rossby waves behave similarly and tend 
to the Alfven-like wave solutions, whereas for n < 3, i.e. for 
large spatial scales, the two modes behave differently: the solu- 
tion with higher negative frequency corresponds to HD Rossby 
waves (dashed line) and the low frequency (continuous line) so- 
lution, which differs from pure Rossby wave dispersion curves, 
is a new wave mode (see similar consideration by Hide 1966 
Acheson & Hide |1973l l. It must be mentioned that the value of 
the toroidal wavenumber, s, does not influence the general be- 
haviour of the modes and that increasing s leads only to a linear 
increase of the frequency, a>. 

The frequency of particular harmonics depends on the 
strength of the toroidal magnetic field. The frequency of the 
s = 1 , n -2 harmonics of fast (dashed) and slow (continuous) 
magnetic Rossby modes vs. the ratio of the Alfven speed to the 
rotation rate, a - vaI2Q.qRq, are plotted in Fig. 3. The fast and 
slow modes have lower frequency than the angular velocity, Qo, 
for a < 0.5 and a < 0.8 respectively. When a goes to zero, 
the fast mode frequency tends to 0.3Qo, while the slow mode 
frequency tends to zero. Therefore, the slow magnetic Rossby 
mode may have very low frequency depending on the ratio of 
the Alfven speed to the angular velocity. 

5. Discussion 

The recently developed "shallow water" magnetohydrodynamic 
approximation (Gilman 2000) has stimulated further study of 
MHD wave modes in this system. While ordinary "shallow wa- 
ter" modes (Poincare, Kelvin, Rossby) have been intensively 
studied in the geophysical context (Pedlosky 119871 Gill |1982| l, 
the inclusion of a horizontal magnetic field enriches the wave 
spectrum. Magneto-gravity and Alfven modes in the "shallow 
water" MHD system have been recently studied by Schecter et 
al. (1200 11 1. However, large-scale modes (those with stellar spatial 
dimensions; for example, Rossby waves) were absent from their 
consideration due to the use of the /-plane approximation. On 
the contrary, here we emphasize the large-scale behavior of wave 
modes in the "shallow water" system. Considering the yS-plane 
approximation and, especially, spherical coordinates, enables us 
to study the wave dynamics on spatial scales corresponding to 
stellar dimensions. Particular attention is paid to the slow mag- 
netic Rossby mode, expressed by Eqs. (20), (45), (48). The pres- 
ence of this mode in rotating magnetised fluids was first pointed 
out by Hide (1966) using the two dimensional Cartesian yS-plane 
approximation in the context of the Earth's liquid core. Here 
we derive analytical dispersion relations of this mode in both 
Cartesian and spherical "shallow water" MHD systems. In the 
low Alfven speed limit (compared to the surface gravity speed), 
this mode has a smaller frequency than that of pure Alfven and 
Rossby modes and consequently may have new interesting con- 
sequences in large-scale stellar dynamics. 

However, this consideration needs some modifications when 
applied to concrete astrophysical situations; for example, to the 
solar tachocline. It is believed that the tachocline is divided 
in two parts: the inner "radiative" layer with a strongly stable 
stratification and the outer "overshoot" layer with a weakly sta- 
ble stratification (Gilman 2000). The subadiabatic stratification 
provides negative buoyancy in both layers, which leads to the 
so-called "reduced gravity", g r (Gilman [2000; Schecter et al. 
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120011 ). Therefore, the developed theory of fast and slow magnetic 
Rossby waves should be modified for tachoclines of solar-like 
stars using the reduced gravity instead of the ordinary one. Then 
the results can be quite different for the radiative and overshoot 
layers due to the significant difference between the reduced grav- 
ity there. Schecter et al. (120011) estimated the reduced gravity 
in the radiative layer as 500-1.5 ■ 10 4 cm-s~ 2 and in the over- 
shoot layer as 0.05-5 cm-s~ 2 . Then the surface gravity speed Co 
is higher than the Alfven speed in the radiative layer, but not in 
the overshoot one, where both speeds may have similar values 
(Schecter et al. 2001 ). Therefore, the results obtained here can 
be easily applied to the radiative part of tachocline, but not to 
the overshoot one. The exception is the dispersion relation (15), 
which can be applied to both parts of the tachocline. 

It must be also mentioned that differential rotation, typical of 
the tachocline dynamics, is absent from our consideration. The 
goal of this paper is to study the influence of the magnetic field 
on the large-scale dynamics of a "shallow water" system in gen- 
eral. However, differential rotation should be taken into account 
in more realistic models of wave dynamics in the tachocline. 

6. Conclusions 

The influence of a toroidal magnetic field on the dynamics of 
Rossby waves in a thin layer of ideal conductive fluid on a ro- 
tating sphere is studied in the shallow water MHD approxima- 
tion for both, Cartesian and spherical geometries. It is shown 
that in both cases the magnetic field causes the splitting of low 
order (long wavelength) ordinary Rossby wave harmonics into 
two modes (here called magnetic Rossby modes). The high fre- 
quency mode (the fast magnetic Rossby mode) corresponds to 
ordinary HD Rossby waves slightly modified by the magnetic 
field, while the low frequency solution leads to a new mode (the 
slow magnetic Rossby mode) with interesting properties. Low 
order (with respect to the poloidal wavenumber) harmonics of 
the slow magnetic Rossby mode have lower frequency than the 
pure Rossby and Alfven wave frequencies of the corresponding 
harmonics. 
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